Dynamic demonstration method and system for water-soluble ion concentration and composition of aerosol

ABSTRACT

Disclosed is a dynamic demonstration method for water-soluble ion concentration and components of an aerosol. The method comprises: obtaining concentration data of each ion in an atmospheric aerosol of a target city in a preset time period and filling the concentration data in a data table; obtaining vertex coordinates of each ion in a Maucha graph according to equivalent concentration data of each ion; drawing an aerosol ion Maucha graph of the target city in each preset time period according to the vertex coordinates; and finally making a dynamic picture according to a temporal graph of aerosol ion concentration in various time periods.

BACKGROUND OF THE INVENTION 1. Technical Field

The present disclosure relates to the field of scientific measurement,and in particular to a dynamic demonstration method and system forwater-soluble ion concentration and composition of an aerosol.

2. Description of Related Art

Aerosol refers to a relatively stable suspension system formed by solidor liquid particles uniformly dispersed in a gas. The clouds in the sky,the morning mist in the forest, the rising sand and dust, the dropletsfrom the ocean, the smoke and dust emitted by industries, the dustformed during construction on construction sites, pesticide spraying,artificial rainfall, etc. are all concrete manifestations of an aerosolin life. Atmospheric aerosol can affect solar radiation and theformation of cloud condensation nuclei and may form photochemicalpollution under light conditions; moreover, fine particles in theaerosol can also enter human lungs through the respiratory tract anddeposit in the alveoli, causing organ diseases. It can be seen that theaerosol has a very important impact on climate change, atmosphericenvironment and human health. As more and more attention is paid to theproblem of air pollution in China, the research on the atmosphericaerosol has gradually become a hot topic in scientific research.

The chemical composition of the atmospheric aerosol is very complex, andwater-soluble inorganic ions are an important part of the atmosphericaerosol. The water-soluble components include a variety of importantsubstances, such as sulfates, nitrates, and the like. These componentscan change the size, composition, pH, quantity and life of the aerosolthrough moisture absorption. In addition, the water-soluble componentscan also increase the solubility of toxic and harmful substances andcause harm to human health. Inorganic ions in the aerosol can be dividedinto primary ions and secondary ions according to their sources. Theprimary ions refer to inorganic ions directly discharged into theatmosphere by an emission source, and the secondary ions refer toinorganic ions converted from the primary ions through chemicalreactions in the atmosphere. Sulfates, nitrates, and ammonium salts,also known as SNA, are important inorganic ions in the secondary aerosoland are important indicators to characterize regional pollution. Theseions mainly come from secondary aerosols formed during thetransformation of gas particles in the atmosphere.

Research on inorganic ions in the aerosol helps to reveal the source andcomposition of the aerosol and is an important method for studying theatmospheric aerosol. At present, monitoring stations have been set upall over the country to monitor the changes in the concentration ofvarious pollutant components in the atmosphere. As a result, a largeamount of monitoring data has been accumulated for people to analyze.However, there is still a lack of effective methods for processing andmining monitoring data. In addition, due to the mobility of the aerosolitself, it may affect the air quality in many areas for a period oftime. By comparing the changes of water-soluble ions in aerosols indifferent regions, it may help us to understand the source and formationof various components in the aerosols. However, there is still a lack ofsimilar methods to solve this problem. In view of the above two points,a method is provided to dynamically display the relative composition ofwater-soluble ions in aerosol over time on a spatial scale. Bydisplaying the changes of water-soluble ion concentrations in aerosolsin different cities on a map, the changes of water-soluble ionconcentration in different cities can be observed within a period oftime, and then some of connections hidden therein can be furtherrevealed.

BRIEF SUMMARY OF THE INVENTION

In order to solve the above problems and fully display the change in ionconcentration of the atmospheric aerosol according to time change andregional change, the present disclosure provides a dynamic demonstrationmethod for water-soluble ion concentration and composition of anaerosol, including the following steps:

S1: obtaining concentration data of each ion in an atmospheric aerosolof a target city in a preset time period, converting the concentrationdata of each ion into equivalent concentration data of each ion, andfilling the equivalent concentration data of each ion in a data table;

S2: importing the data table, and obtaining a radius of a circle in aMaucha graph and a diagonal length corresponding to a quadrilateral ofeach ion in the Maucha graph according to a first calculation formula onthe basis of the equivalent concentration data of each ion;

S3: obtaining vertex coordinates of each ion in the Maucha graphaccording to a second calculation formula on the basis of the diagonallength;

S4: drawing an aerosol ion Maucha graph of the target city in eachpreset time periods according to the vertex coordinates of each ion inthe Maucha graph; and

S5: determining whether the aerosol ion Maucha graphs of all targetcities are completed or not, if not, switching the target city andreturning to step S1;

wherein the first calculation formula is:R ²×sin22.5°/2=T/16, b×R×sin22.5°/2=P/2;

the second calculation formula is:bx=b×cos(22.5°+45°×n), by=b×sin(22.5°+45°×n);

in the formulas, R refers to the radius of the circle in the Mauchagraph; T refers to total equivalent concentration of eight major ions inthe atmospheric aerosol; b refers to a diagonal length corresponding toa quadrilateral of a corresponding ion in the Maucha graph; P refers tothe equivalent concentration of the corresponding ion; bx refers to theabscissa of a vertex of the corresponding ion; by refers to an ordinateof the vertex of the corresponding ion; n refers to a constant variablethat changes with an angle of a line of hexadec-section in thequadrilateral corresponding to the ion in the circle, and the value of nbegins to vary counterclockwise with the starting angle (22.5°) of theline of hexadec-section and increases from 0 to 7.

Further, the data table also includes time information and latitude andlongitude information of the target city corresponding to theconcentration data of each ion.

Further, before step S1, the method further includes:

Step S0: converting ion categories in a water body in an original Mauchagraph to ion categories in the atmospheric aerosol;

the ion categories in the atmospheric aerosol include eight major ionswhich are K⁺, Na⁺, Ca²⁺, NH₄ ⁺, SO₄ ²⁻, Cl⁻, NO₃ ⁻, and F⁻,respectively.

Further, after step S5, the method further includes:

S6: superimposing the aerosol ion Maucha graph of each target city inthe same time period on a geographical map according to the latitude andlongitude information, and drawing a temporal graph of aerosol ionconcentration; and

S7: making a dynamic picture according to the temporal graph of aerosolion concentration in each time period.

Further, R-Shiny is used for writing and packaging in steps S2 to S7 andpackaging results are displayed on a web terminal.

Further, in step S2, in the process of importing the data table, thedata need to be determined. If there are more than a preset number ofmissing data in a row of the data table, the row of data is deleted.

The present disclosure further provides a dynamic demonstration systemfor water-soluble ion concentration and composition of an aerosol,including an import module, a calculation module, and a drawing module,wherein

the import module is configured to import equivalent concentration dataof each ion in an atmospheric aerosol of a target city in a preset timeperiod in a data table and transmit the equivalent concentration data tothe calculation module;

the calculation module is configured to obtain vertex coordinates ofeach ion in a Maucha graph according to a first calculation formula anda second calculation formula on the basis of the equivalentconcentration data of each ion;

the drawing module is configured to draw an aerosol ion Maucha graph ofthe target city in each preset time period according to the vertexcoordinates of each ion in the Maucha graph;

wherein the first calculation formula is:R ²×sin22.5°/2=T/16, b×R×sin22.5°/2=P/2;

the second calculation formula is:bx=b×cos(22.5°+45°×n), by=b×sin(22.5°+45°×n);

in the formulas, R refers to the radius of the circle in the Mauchagraph; T refers to the total equivalent concentration of eight majorions in the atmospheric aerosol; b refers to a diagonal lengthcorresponding to a quadrilateral of the corresponding ion in the Mauchagraph; P refers to the equivalent concentration of the correspondingion; bx refers to the abscissa of a vertex of the corresponding ion; byrefers to an ordinate of the vertex of the corresponding ion; n refersto a constant variable that changes with an angle of a line ofhexadec-section in the quadrilateral corresponding to the ion in thecircle, and the value of n begins to vary counterclockwise with thestarting angle (22.5°) of the line of hexadec-section and increases from0 to 7.

Further, the data table also includes time information and latitude andlongitude information of the target city corresponding to theconcentration data of each ion.

Further, the import module is further configured to convert ioncategories in a water body in an original Maucha graph to ion categoriesin the atmospheric aerosol;

the ion categories in the atmospheric aerosol include eight major ionswhich are K⁺, Na⁺, Ca²⁺, NH₄ ⁺, SO₄ ²⁻, Cl⁻, NO₃ ⁻, and F⁻,respectively.

Further, the drawing module further includes a superimposing unit,

configured to superimpose the aerosol ion Maucha graph of each targetcity in the same time period on a geographical map according to thelatitude and longitude information, draw a temporal graph of aerosol ionconcentration, and make a dynamic picture according to the temporalgraph of aerosol ion concentration in each time period.

Further, the import module further includes a determining unit,

configured to determine the data, wherein if there are more than apreset number of missing data in a row of the data table, the row ofdata is deleted.

Further, R-Shiny is used for writing and packaging in the import module,the calculation module, the drawing module, the superimposing unit, andthe determining unit, and packaging results are displayed on a webterminal.

Compared with the prior art, the present disclosure at least has thefollowing beneficial effects:

(1) The dynamic demonstration method and system for water-soluble ionconcentration and composition of an aerosol in the present disclosuredisplays the water-soluble ions in the atmospheric aerosol through theMaucha graph, so that the concentration of each ion can be displayedmore intuitively.

(2) By integrating the aerosol ion Maucha graphs of various cities inthe same time period on the same geographical map, the temporal graph ofaerosol ion concentration is made, and the temporal graphs in varioustime periods are integrated into a dynamic picture, which can moreintuitively display the concentration change trend and composition ofwater-soluble ions in the atmospheric aerosol of different cities indifferent time periods in the same time dimension.

(3) R language (Shiny data package) is used to package and writeprograms, and the corresponding tasks (such as data import, dataanalysis, drawing, etc.) can be completed according to simple datawithout relying on front-end and back-end engineers. Because the Rlanguage is used and the operation can be easily displayed on the webterminal, it is easier to be accessed by other users, thus facilitatingthe transmission of information.

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS

FIG. 1 is a schematic diagram of method steps of a dynamic demonstrationmethod for water-soluble ion concentration and composition of anaerosol;

FIG. 2 is a schematic diagram of a system structure of a dynamicdemonstration system for water-soluble ion concentration and compositionof an aerosol;

FIG. 3 is a schematic diagram of composition of a Maucha graph; and

FIG. 4 is a schematic diagram of a Maucha graph drawing method.

DETAILED DESCRIPTION OF THE INVENTION

The following are specific embodiments of the present disclosure andfurther describe the technical solutions of the present disclosure inconjunction with the accompanying drawings, but the present disclosureis not limited to these embodiments.

Embodiment 1

In order to solve the above problems and fully display the change in ionconcentration of the atmospheric aerosol according to time change andregional change, as shown in FIG. 1 , the present disclosure provides adynamic demonstration method for water-soluble ion concentration andcomposition of an aerosol, including the following steps.

In step S1, concentration data of each ion in an atmospheric aerosol ofa target city in a preset time period is obtained, the concentrationdata of each ion is converted into equivalent concentration data of eachion, and the equivalent concentration data of each ion is filled in adata table.

In the data table, in order to obtain more comprehensive information forsubsequent drawing, the data table further includes corresponding timeinformation, longitude and latitude information of the target citycorresponding to the concentration data of each ion, and the like, sothat areas can be divided according to the longitude and latitude duringthe subsequent drawing.

Certainly, the Maucha graph should not be simply used directly. Becausethe original Maucha graph is mainly used for water bodies, but in thepresent invention the original Maucha graph is used for an atmosphericaerosol, the water body and the atmospheric aerosol are different in ioncategories, and the ion categories of the atmospheric aerosol need to beadjusted. Therefore, before step S1, the method further includes:

Step S0: converting ion categories in a water body in an original Mauchagraph to ion categories in the atmospheric aerosol;

the ion categories in the atmospheric aerosol include eight major ionswhich are K⁺, Na⁺, Ca²⁺, NH₄ ⁺, SO₄ ²⁻, Cl⁻, NO₃ ⁻, and F⁻,respectively.

Because an atmospheric combination of the above-mentioned ions containsa variety of chemical pollutants, the concentrations of these ions canbe detected to analyze the concentrations of chemical pollutants such ashydroxides, hydrocarbons, and sulfides in the atmosphere.

In step S2, the data table is imported, and a radius of a circle in aMaucha graph and a diagonal length corresponding to a quadrilateral ofeach ion in the Maucha graph are obtained according to a firstcalculation formula on the basis of the equivalent concentration data ofeach ion.

In some cases, data missing may be detected. If the number of missingdata reaches a preset number (specifically a value that makes a drawingfunction wrong, for example, it is assumed that the preset number is setto 2, there are 7 data units when all data in a row is not missing; whenin the imported data table, there are missing data in 2 or more dataunits in a row, then the drawing function will be biased), in order toavoid the adverse effect of missing data on data analysis and drawing,the process of importing the data table in step S2 further includes datadetermination. If there are more than a preset number of missing data ina row of the data table, the row of data is deleted.

As shown in FIG. 3 , the Maucha graph consists of a circle and 8quadrilaterals each having two sides located in the circle. The area ofeach quadrilateral represents the equivalent concentration of an ioncorresponding to the quadrilateral in the atmosphere. Specifically, thecircle is equally divided into four ¼ circles up, down, left, and right.Each ¼ circle is divided into two ⅛ circles. Eight quadrilaterals arelocated within the eight ⅛ circles and their extension areas. Two sidesof the ⅛ circle connected to the center of the circle constitute twosides of the quadrilateral, and the center of the circle serves as oneof the corner points of each quadrilateral. The radius of the circle isdetermined by the area of an inscribed regular hexadecagon of thecircle, and the area of the inscribed regular hexadecagon should beconsistent with the total concentration of the eight major ioncomponents mentioned above. The quadrilateral representing a certain ioncomponent is determined as follows. As shown in FIG. 4 (in this figure,Point O represents the center of the circle at point (0,0) on an xycoordinate axis; Points D, C G, H, I, J, K, and L are oct-section pointsof the circle, lines of oct-section are line segments that connect thecenter of the circle and the oct-section points of the circle, that is,the line segments OD, OC, OG, OH, OI, OJ, OK, OL, the lines ofoct-section divide the circle into eight equal parts; Point A representsa hexadec-section point of the circle on a CD Arc, lines ofhexadec-section are line segments that connect the center of the circleand the hexadec-section points of the circle, the lines ofhexadec-section divide the circle into eight sixteen equal parts; PointB is hypothetical and represents a vertex of a certain ion on a line ofhexadec-section at a certain equivalent concentration of the ion; PointE represents a vertical point of Point A on a line segment OC; Point Frepresents a vertical point of Point C on a line segment OA): when thecircle is equally divided into eight sectors, for each sector, aquadrilateral can be defined by two line segments (such as the linesegment OD and the line segment OC) that make up the sector and a point(e.g., Point B, which is on a ray passing a hexadec-section point of thecircle) inside or outside the sector, and the area of the quadrilateral(shaded part as shown in FIG. 3 ) is required to be equal to theequivalent concentration of the ion. In this way, another 7quadrilaterals representing other ions are sequentially determined, andthe sum of the areas of these eight quadrilaterals must be equal to thearea of the inscribed regular hexadecagon of the circle. Therefore, thevertex coordinates of each ion can be obtained by determining the radiusof the circle and the corresponding ion concentration to draw a Mauchagraph.

In step S3, vertex coordinates of each ion in the Maucha graph areobtained according to a second calculation formula on the basis of thediagonal length.

In step S4, an aerosol ion Maucha graph of a target city in each presettime period is drawn according to the vertex coordinates of each ion inthe Maucha graph.

In step S5, it is determined whether the aerosol ion Maucha graphs ofall target cities are completed or not; if not, the target city isswitched and the process returns to step S1.

wherein the first calculation formula is:R ²×sin22.5°/2=T/16, b×R×sin22.5°/2=P/2;

the second calculation formula is:bx=b×cos(22.5°+45°×n), by=b×sin(22.5°+45°×n);

In the formula, R refers to the radius of the circle in the Mauchagraph; T refers to the total equivalent concentration of eight majorions in the atmospheric aerosol; b refers to the diagonal lengthcorresponding to the quadrilateral of the corresponding ion in theMaucha graph (as shown in FIG. 3 and FIG. 4 , the diagonal correspondingto the quadrilateral of the corresponding ion in the Maucha graph isdefined as the line segment connecting the center of the circle—Point Oand the vertex—e.g., Point B); P refers to the equivalent concentrationof the corresponding ion; bx refers to the abscissa of the vertex of thecorresponding ion; by refers to the ordinate of the vertex of thecorresponding ion; n refers to a constant variable that changes with anangle of a line of hexadec-section in the quadrilateral corresponding tothe ion in the circle, and the value of n begins to varycounterclockwise with the starting angle (22.5°) of the line ofhexadec-section and increases from 0 to 7. For example (as shown in FIG.4 ), when the quadrilateral corresponding to the current ion is in anarea defined by COG, that is, the line of hexadec-section in thequadrilateral is the starting angle of 22.5°, n is 0, and when thequadrilateral corresponding to the current ion is in the area defined byCOG, that is, the line of hexadec-section in the quadrilateral is theangle of 45°, n is 1.

In the meanwhile, after step S5, in order to make the final image to beable to fully show the change and transfer trend of ion concentration inthe atmospheric aerosol, the following steps are also included.

In step S6, the aerosol ion Maucha graph of each target city in the sametime period is superimposed on a geographical map according to thelatitude and longitude information, and a temporal graph of aerosol ionconcentration is drawn.

S7: making a dynamic picture according to the temporal graph of aerosolion concentration in each time period.

By superimposing the aerosol ion concentration Maucha graph of each cityat the same time on the geographical map, the distribution of pollutantsin various places can be clearly shown, and a dynamic picture is madetherefrom. On the original basis, the changes of pollutants in variousplaces can be reflected, and the diffusion and movement of pollutantscan be analyzed according to the increase or decrease of ionconcentration in various places, which is conducive to the tracking ofpollutant sources.

It should be noted that R-Shiny is used for writing and packaging insteps S2 to S7 and packaging results are displayed on a web terminal.R-Shiny is used for writing and packaging because Shiny is a webdevelopment framework of R language by which users only need tounderstand some html knowledge to quickly complete web developmentwithout deep understanding of css and js. Moreover, the shiny packageintegrates features such as bootstrap, jquery, ajax, etc., which greatlyliberates the productivity of R as a statistical language. In this way,R users of non-traditional programmers can complete some simple datavisualization tasks according to their business without relying onfront-end and back-end engineers, and quickly verify the reliability oftheir ideas.

The Shiny application contains two basic components: one is a userinterface script, and the other is a server script. In this method, astreaming page in a user interface is mainly used to put parts such asdata processing, drawing, and result presentation into the serverscript. A series of operations such as data import, data analysis, andgraph drawing are packaged in a shiny package to form a web application,so that the web application can be opened in a browser and the finalresults are presented on the network.

According to the generalized overview, the data table is importedthrough Shiny (Shiny is an application framework that provides a webinteractive interface for an R model) and whether data in the data tableis complete is determined; then, the data are analyzed to draw graphsand the obtained aerosol ion Maucha graph of each city in each timeperiod is saved in a first folder (set as Folder figures) through a pngdata package; the aerosol ion Maucha graph of each city in the same timeperiod is superimposed on a corresponding position of the geographicalmap (the geographical map is a satellite map, provided by anotherfolder, set as Folder maps) through a ggimage data package, thusobtaining the temporal graph of aerosol ion concentration of each cityin a single time period on the geographical map; and the temporal graphis saved in a second folder (set as Folder figures2) Then, when the uservisits a web segment, the temporal graph of the aerosol ionconcentration in the corresponding time period can be called through amagick data package and made into a gif dynamic picture, and dataanalysis and trend determination are carried out on the basis of thepicture. It should be noted that the calling of the above-mentioned datapackets can be implemented through a getGIF.r script file, and finally aroute of the gif picture is returned to the shinyAPP server script, andthe shinyAPP server script calls the picture according to the route.

Embodiment 2

In order to better analyze the inventive idea of the present disclosure,this embodiment divides the functions described in the invention in aform of system structure. As shown in FIG. 2 is a dynamic demonstrationsystem for water-soluble ion concentration and composition of anaerosol, including an import module, a calculation module, and a drawingmodule, wherein

the import module is configured to import equivalent concentration dataof each ion in an atmospheric aerosol of a target city in a preset timeperiod in a data table and transmit the equivalent concentration data tothe calculation module;

the calculation module is configured to obtain vertex coordinates ofeach ion in a Maucha graph according to a first calculation formula anda second calculation formula on the basis of the equivalentconcentration data of each ion;

the drawing module is configured to draw an aerosol ion Maucha graph ofthe target city in each preset time period according to the vertexcoordinates of each ion in the Maucha graph;

wherein the first calculation formula is:R ²×sin22.5°/2=T/16, b×R×sin22.5°/2=P/2;

the second calculation formula is:bx=b×cos(22.5°+45°×n), by=b×sin(22.5°+45°×n);

in the formulas, R refers to the radius of the circle in the Mauchagraph; T refers to the total equivalent concentration of eight majorions in the atmospheric aerosol; b refers to a diagonal lengthcorresponding to a quadrilateral of the corresponding ion in the Mauchagraph; P refers to the equivalent concentration of the correspondingion; bx refers to an abscissa of a vertex of the corresponding ion; byrefers to an ordinate of the vertex of the corresponding ion; n refersto a constant variable that changes with an angle of a line ofhexadec-section in the quadrilateral corresponding to the ion in thecircle, and the value of n begins to vary counterclockwise with thestarting angle (22.5°) of the line of hexadec-section and increases from0 to 7.

The data table also includes time information and latitude and longitudeinformation of the target city corresponding to the concentration dataof each ion, and the corresponding data is input by an input module(usually input tools such as a keyboard).

The import module is further configured to convert ion categories in awater body in an original Maucha graph to ion categories in theatmospheric aerosol;

the ion categories in the atmospheric aerosol include eight major ionswhich are K⁺, Na⁺, Ca²⁺, NH₄ ⁺, SO₄ ²⁻, Cl⁻, NO₃ ⁻, and F⁻,respectively.

The drawing module further includes a superimposing unit,

configured to superimpose the aerosol ion Maucha graph of each targetcity in the same time period on a geographical map according to thelatitude and longitude information, draw a temporal graph of aerosol ionconcentration, and make a dynamic picture according to the temporalgraph of aerosol ion concentration in each time period.

The import module further includes a determining unit,

configured to determine the data, wherein if there are more than apreset number of missing data in a row of the data table, the row ofdata is deleted.

R-Shiny is used for writing and packaging in the import module, thecalculation module, the drawing module, the superimposing unit, andpackaging results are displayed on a web terminal (through a displaymodule, i.e., a display).

In summary, the dynamic demonstration method and system forwater-soluble ion concentration and composition of an aerosol in thepresent disclosure displays the water-soluble ions in the atmosphericaerosol through the Maucha graph, so that the concentration of each ioncan be displayed more intuitively. By integrating the aerosol ion Mauchagraphs of various cities in the same time period on the samegeographical map, the temporal graph of aerosol ion concentration ismade, and the temporal graphs in various time periods are integratedinto a dynamic picture, which can more intuitively display theconcentration change trend and composition of water-soluble ions in theatmospheric aerosol of different cities in different time periods in thesame time dimension.

R language (Shiny data package) is used to package and write programs,and the corresponding tasks (such as data import, data analysis,drawing, etc.) can be completed according to simple data without relyingon front-end and back-end engineers. Because the R language is used andthe operation can be easily displayed on the web terminal, it is easierto be accessed by other users, thus facilitating the transmission ofinformation.

The specific embodiments described herein are only examples toillustrate the spirit of the present disclosure. Those skilled in theart to which the present disclosure belongs can make variousmodifications or additions to the specific embodiments described or usesimilar alternatives, without departing from the spirit of the presentdisclosure or going beyond the scope defined in the appended claims.

What is claimed is:
 1. A dynamic demonstration method for ionconcentration and composition of an aerosol, comprising the followingsteps: S0: labeling quadrants of a Maucha graph with eight major ions inan atmospheric aerosol, wherein the eight major ions in the atmosphericaerosol are K⁺, Na⁺, Ca²⁺, NH₄ ^(+, SO) ₄ ²⁻, Cl⁻, NO₃ ⁻, and F⁻,respectively; S1: obtaining concentration data of each ion in theatmospheric aerosol of one of a plurality of target cities in eachpreset time period, converting the concentration data of each ion intoequivalent concentration data of each ion, and filling the equivalentconcentration data of each ion in a data table; S2: obtaining a radiusof a circle in a Maucha graph according to a first calculation formula,and obtaining a diagonal length corresponding to a quadrilateral of eachion in the Maucha graph according to the first calculation formula onthe basis of the equivalent concentration data of each ion in the datatable and the radius of the circle in the Maucha graph; S3: obtainingvertex coordinates of each ion in the Maucha graph according to a secondcalculation formula on the basis of the diagonal length; S4: drawing anaerosol ion Maucha graph of the one of a plurality of target cities ineach preset time period according to the radius of a circle in a Mauchagraph and the vertex coordinates of each ion in the Maucha graph; andS5: determining whether the aerosol ion Maucha graphs of all of aplurality of target cities are completed or not, and if not, switchingthe target city and returning to step S1; wherein the first calculationformula is:R ²×sin22.5°/2=T/16 b×R×sin22.5°/2=P/2; the second calculation formulais:bx=b×cos(22.5°+45°×n), by=b×sin(22.5°+45°×n); in the formulas, R refersto the radius of the circle in the Maucha graph; T refers to totalequivalent concentration of the eight major ions in the atmosphericaerosol; b refers to a diagonal length corresponding to a quadrilateralof a corresponding ion in the Maucha graph; P refers to the equivalentconcentration of the corresponding ion; bx refers to the abscissa of avertex of the corresponding ion; by refers to an ordinate of the vertexof the corresponding ion; n refers to an integer from 0 to 7; S6:superimposing the aerosol ion Maucha graphs of each of the plurality oftarget cities in the same preset time period on a geographical mapaccording to the latitude and longitude information, to form a temporalgraph of aerosol ion concentration; and S7: repeating step S6 for eachof the preset time periods; S8: superimposing the temporal graphs ofeach aerosol ion concentration in each preset time period to form adynamic picture.
 2. The dynamic demonstration method for ionconcentration and composition of an aerosol according to claim 1,wherein the data table also comprises time information and latitude andlongitude information of the target city corresponding to theconcentration data of each ion.